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ABSTRACT 



SPACE /vND ENERGY DEPENDENT NEUTRON 
SLO’vVIRG DO'.VR TIME IN GRAPHITE 



Submitted to the Departtient of Nuclear Engineering on 
September 26 , 196 ? in partial fulfillment' of the requirements 
for the degree of Master of Science. 

A 4 foot cube, graphit§ moderating assembly was con- 
structed at the M.I.T. Rockefeller Van de Graaff accelerator. 
Using the pulsed neutron experimental installation, a neu- 
tron burst of 200 nsec width was injected into the center ■ 
^of the assembly. Prompt gamma radiation from a neutron 
resonance capture by an indicator foil in the moderator 
was measured, with a 5 inch x 5 inch sodium iodide scin- . 

tillation detector, as a function of time after the neu- j 
tron pulse. The energy dependence of slowing down time { 

was experimentally raeasured by using resonance absorber ; 

foils of cadmium, indium, and gold with isolated reson- 
ances of 0 . 2 , 1.46 and 5*0 eV respectively. By vary- 
ing the distance of the foil from the source, the spatial ' 
dependence of slowing down time was also investigated. 
Comparisons are made with slowing dovm theory. The experi- : 
mental results, together v/ith their estimated uncertainties,' 
are tabulated below: 
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CHAPTER 1 
Introduction 

1 

In his M.I.T. Ph.D. thesis, Adamantiades (l) demon- 

•V 

strated that the pulsed neutron technique using a resonance 
absorber as a point indicator could be used to measure sys- 
tem transient effects such' as the sleeving down time. He 
designed- and installed an experimental installation for 
pulsed .neutron experiments, and made some measurements in 
water. Vaughan (2) refined Adamantaides ’ system, and measured 
the spatial dependence of slowing down time in v/ater for 

three different resonance absorbers. Moller has carried 

u 

out extensive experiments in the field of neutron slowing 
do^vn for both light (5,^) and heavy (5) water. Hov/ever, 
there does not appear to be any published measurements of 
the slowing down time for heavier moderators. Therefore, 
in this present experiment, the pulsed technique has been 
expanded to investigate slov/ing down time in graphite. 

In the pulsed neutron method, a short monoenergetic ' 
burst of neutrons is injected into the moderating medium 
at a point. The neutrons, originally high energy, gradually 
slow down by elastic collisions with the moderating medium. 

If a neutron indicator is incorporated into the moderating 
medium, it is possible to measure the time difference betv/een 
pulse injection and neutron arrival at the indicator. A 
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time distribution is formed from many time differences, 

measured under the exact same conditions. 

! 

Since most materials emit prompt gamma rays upon the 
capture of a neutron, detection of these gamma rays pro- 
(Vides one method by v/hich the time behavior of neutron 
captures may be followed after the pulse. If an absorber 

I 

with a strong, isolated neutron capture resonance is incor- 
porated into the moderating system as a point indicator, 
then the- reaction rate of the neutron flux v/ith the absorp- 
tion resonance is proportional to the time dependent flux 
at that point, for neutrons with energies equal to the 
resonance energy. That is, the neutron density (or flux) 

N(r,E t) 

iJ . 

N - neutron density 
r - distance from source 

- energy of absorber resonance 
t - time 

can be measured. Thus, variation of the position of the 
absorber and the use of different resonance absorbers should 
yield the time, space, and energy variation of the neutron 
density following the pulsing of the assembly. 

By using the M.I.T. Rockefeller Van de Graaff accelera- 
tor as a source of high energy protons, high energy neutrons 
were obtained by the Li*^(p ,n)Be*^^ nuclear reaction. The lithiinn 
target was located at the center of a large assembly of 
graphite, thus approximating a point source in an infi- 
nite medium. The accelerator beam was pulsed by an external. 



pulsing sysbem after the protons had been accelerated. The 
pulse width used in these experiments v/as 200 nsec. Enough 
time was allowed' to elapse between pulses so that all tran- 
sient neutron effects from one pulse would die out before 
] • 

; * 

the next pulse was injected. Absorber foils of cadmium, 
indium, and gold, with isolated resonances located at 0.2 eV, 
1.46 eV, and 5*0 eV respectively, were used to sample the 
neutron flux at these energies. A detecting and analyzing 
system was assembled that was capable of indicating the 

* 

time difference between pulse injection and the resonance 
capture of a neutron. The most probable time difference 
from injection to resonance capture was assumed to be the 
neutron slowing down time for that energy at that position. 

The material in this thesis is organized as' follows: 
Chapter 2 is concerned with slowing dovm theory^ the , 

experimental facilities are described in Chapter 3 ] the 
performance of the experiments is described in Chapter' 4^ 
and the results are presented and discussed. 



I 
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' . . CHAPTER 2 ' 

Theory 

’ 

2,1 General ^ 

The analytical treatment of pulsed neutron assem- 
■ hlies requires a mathematical solution of the general 
transport equation. If an' expression representing the neu- 
tron density as a function of time, energy, and space, and 
satisfying the Boltzmann transport equation with its source 
and boundary conditions can be found, the behavior of neu- 
trons in the 'system can be completely described. 

.• t 

Unfortunately, no exact solution of the Boltzmann 

! j 

equation, in all its generality, has ever been found^ 
but .approximate solutions are available. Of these solutions, 
, only-a few may be verified because not all quantities can 
be.'measured experimentally. 

The quantity most readily measured experimentally 
in the slowing down region is the neutron slowing down 
time (tj^) . This is the time at v/hich the neutron density 
at some energy (E) is maximum and can be found by setting 
the derivative of N(r,E,t) 

N - neutron density 
r - distance from source 

• % 

E - neutron energy 



\ 



t - time 
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with respect to time equal to zero. 

It is also possible to measure the mean slowing dov/n 
time (<t>) which is also defined at the first time moment. 
The tim*e moments are integral quantities defined as: 

<t^(r,v)> N(v,t) dt/f‘^N(v,t) dt (2.1.1) 

lT(v,t) - neutron density as a function of space 
and energy (velocity) 
t - time 

Both the time to maximum or the slov/ing down time, 
and mean slov/ing dov/n time may be measured experimentally as 
a function of space and neutron energy. 

In this experiment, the source may be considered 
isotropic and a delta function in time, energy, and space. 
Although the neutron source energy does vary v/ith the angle 
in the laboratory system, this source term can be neglected' 
because it appears in the denominator of the expressions 
for the time distributions, and it is much greater than the 
energy at which the time distribution is measured. For 
example, one expression for the neutron slowing down flux 
in hydrogen is: 

(p(v,t) = S.(vtAg)" *exp(-vt/Yg) *(l-v/v^+2Yg/v^t) + 
2S*v^/v-6(v-v^) -exp(-v^t/Yg) (2.1.2) 
cp - neutron flux 
.V - neutron velocity 
S - source strength 



Of 
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Y„ - scattering mean free path 

O ' * 

- neutron source velocity 
t - time 

However, because v^»v, equation (2.1.2) simplifies to: 
cp(v,t) =■ S- (Jt/Yg)^ * exp(-vt/Yg) (2.1.5) 

Except for positions closer than three mean free 
paths away from the physical boundary of the assembly, the 
spatial boundary conditions are those of an infinite medium. 

9 

For positions closer than three mean paths, leakage must 
be taken into account. 

In the energy region lying above a few tenths of an 
eV, the nuclei of the moderating assembly can be considered 
free and at rest. For lower energies, hov/ever, the molecular 

It 

structure of the moderator becomes important. 

Any absorption v/hich occurs in the moderating assembly 
will effect the time distribution. Those neutrons which 
take longer to slow down, i.e. have more collisions, to 
energy (v) v/ill have a greater probability of being absorbed. 
Thus, absorption contributes to a decrease in t . In 
graphite, the absorption cross section is very small^ and 
absorption effects may be neglected. 

2.2 Space Independent Slowing; Down Theory 

If, in addition to an infinite medium, we assume a 
uniformly distributed source instead of a spatial delta 
function source, the spatial dependence is removed from the 
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Boltzm^mn Equation. The equation reduces to: 

. 1/v- 6(p(v, t)/6t = -[e (v) + e (v)](p(v,t) + 

o ci 

■ Jcp(v’ , t) • e( v’ — >v) dv’+S*6(v-v^) ‘SCt) (2.2.1) 

(p(v,t) - neutron flux as a function of space and 
energy (velocity) 

• , e - macroscopic scattering cross section 
1 s 

- macroscopic absorption cross section 
e - removal cross section 
s -.source strength (n/cm*-sec) 
v^ - source energy 
6( ) - Dirac delta function 

Unfortunately, for a heavier than water moderator, 

it is very difficult to obtain an expression for the flux, 

i) 

and hence, for the slowing down time. Using a method 
introduced by Marshak (1) the time moments can be obtained 
exactly. A few of the moments, in terms of the moderator 
mass number, obtained by Adamantiades (2) are given in 
Table I, of his thesis. Once the time moments have been 
determined, trial functions can be found to fit these 
moments . 

Marshak and vonDardel (3) suggest using the function: 
cp(x) = ^ • exp(-x-b/x) • (2.2.2) 

A - normalization constant 

X = vt/y where y - scattering mean free path 
r = (M-m)/(M+m) 

M - mass of moderator 



m - mass of neutron 
b - a parameter 

which gives the right asymptotic dependence for the higher 
positive moments. vonDardel varies the parameter (b) with 

S 

(r) while Marshak treats it as a constant. 

■ I 

j. For heavier elements, the time dependent neutron den- 

sity distribution approaches a narrow Gaussian in shape. 
Thus, for heavier elements, t^ should approximately equal 
<t>. 

Adamantiades gives a complete discussion of the time 
moment method in Chapter II of his thesis. He also men- 
tions other methods, and gives a number of references deal- 
ing with the problem of a heavier than water moderator. 

l; 

None of these references, ho'wever, specifically give result 
for graphite assemblies. 

2.3 Space Dependent Slowing; Dovm Theory 

When the neutron source is localized, as in this 
experiment, the effects of diffusion during moderation must 
be taken into account. The problem becomes much more com- 
plicated since the transport equation is now a function of 
three, rather than two variables. 

As in the space independent case, Adamantiades in Chap 
ter II of his thesis discusses the theoretical v/ork which 
has been done in this area. This concludes a solution by 
Dyad ’.kin and Batalina (A), for the space-energy-time neu- 
tron density distribution assuming no absorption, constant 





i 



scattering cross section, and isotropic scattering. The 
expression involves three factors which are complicated^ 
and require coefficients which are dependent on the medium 
being examined. The authors also give a simplified expres- 
sion for t , which is: 
m 

^max " 2Yg/vTi- (2-i3p/YFTFir) . . • 

[l-2a/3*(l-pp/Vr^+r p' (2.3.1) 

I ■ ■ 

Y - scattering mean free path 
s 

V - neutron velocity , . 

Tt - average logarithmic energy decrement 

p - neutron lethargy 

a - parameter which depends on the moderating medium 

P - parameter which depends on the moderating medium 

ii 

r - distance from source 
Results are given for sand and water mixtures. 

Claesson (5) has also determined the space dependent 
neutron density distribution for both a point and plane 
source assuming zero absorption, isotropic scattering, 
and constant mean free path. Although his method can be 
extended to moderators of arbitrary mass, curves are only 
given for the case of water at the indium resonance. These 
curves show however, that past 10 cm the slowing down time 
is strongly dependent on neutron energy at birth. 

Adamantiades also develops a method of calculating 
time moments using an iterative solution of the stationary 
transport equation. The nth time moment is expressed 
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as : 

. <t’"> = ,(2.5.2) 

= -S(r,v,VL) (2.3.5) 

Hcp^^^+n/v* = 0 (2.3.^) 

H = [-y,*'^-e+.+ydv’>/dvL ’ e^(v’ — >v,y,’ — >y.) ] (2.5.5) 

b S 

. r,v,]i - vectors for position, velocity and solid 
\ 

angle, respectively 

- total macroscopic absorption cross section. 

Any calculational method suitable to solve the steady state 
can be used to solve for the flux moments. 

David Diamond, a graduate student at M.I.T., has 
kindly supplied -theoretical time distributions at the gold 

and indium resonances for a number of spatial points, 

i) 

using the diffusion approximation. A distribution for each 
foil at a distance from the source of 18 cm is shown in 
Figure 2.1 (a) and (b) . 



- 16 - 




r-( 






q.qSis;,^ OAi^^Bieg 



Figure 2.1 (a) Theoretical Time Distribution at 1.46 eV 
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CHAPTER 3 

Experimental Equipment 



3.1 General Description 

The equipment necessary for measuring neutron transient 
effects by the pulsed source method can be grouped into six 
categories; 

1) The neutron source. 

‘ ^ 2) A method for pulsing the source. 

' 5) The moderating assembly. 

4) A system for -detecting, analyzing and storing data. 

5) A system to monitor the experimental equipment. 

6) A normalizing system to insure consistent data. 

The general layout ■ of the experimental system showing •' 

I 

the major components is given in Figure 3.1* 

The neutron source consisted of a target of natural 
lithium, located approximately in the center of the modera- 
ting assembly, bombarded by 2.4 MeV protons to produce the 
Li'(p,n)Be'^ reaction. The high energy protons were produced 
by’ a Van de Graaff electrostatic accelerator. The pulsing 
system electrostatically deflected the proton 'beam off the 
lithium target, allov/ing it to strike the target for a short 
time. 

Thus, the pulse of fast neutrons produced at the target 
was ejected into the center of the moderating assembly. 
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Figure 3,1 Block Diagram of Major Systems 
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These rieutrono were then slowed down by elastic collisions 
in the medium, and their time distribut *,on was measured 
by the detecting and analyzing system. This system con- 
sisted of the conventional detector-preamplifier-amplifier 
array, plus electronic equipment for pulse height discri- 
mination and time analysis combined with a multi-channel 
ainalyzer to store the data. There was a common link between 

the pulsing and the analyzing system which provided a zero 

" / 
time reference. _ j 

The monitoring system was used to control the perfor- / 

\ 

mance of the pulsing system, and thus to optimize the quality 
of the neutron pulse. The normalizing system yielded the j 
time integrated neutron flux for each experiment performed. • 



This information was necessary to properly carry out back- 



ground subtraction. 

A more detailed description of the experimental equip- 
ment will be given in the sections to follow. 

I 

3.2 The Mechanical Layout 

A plan view of the mechanical layout of the system is 
given in Figure 3»2. Only a brief description of the sys- 
tem will be given here since a thorough, detailed descrip- 
tion of the system v/ill be given by Adamantiades (1), who 
actually constructed the system. 

Since the experimental apparatus, v/as positioned on 
the -floor below the accelerator, the protons were accelerated 
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dov/n-.'/ard and bent into the horizontal beam tube by a bend- 
ing and analyzing magnet. The beam tube was Icept under a 
high vacuum by the diffusion pump of the accelerator, and 
by a small Vac Ion pump located close to the target. 

The deflecting plates delivered a high voltage pulse 
v/hich deflected the beam. The water-cooled chopping slits 
were mounted after the deflection plates for the purpose 
of intercepting the beam when it was sv/ept off the target. 
The distance betv/een the slits, the shape of the high vol- 
tage pulse, and the deflecting plates configuration deter- 
mined the shape of the pulse reaching the target. 

A remotely controlled pneumatic shutter valve was 
interposed between the chopping slits and the target. 

ii 

The beam could thus be cut off from the target, enabling 
v;ork to be done around the assembly. The target tube itself 
v/as constructed of aluminum and was 23 inches long, 3/8 
inch O.D., with a 1/16 inch wall thickness. The target 

i 

assembly v/as located at the end of the target tube. 

It was decided to cool the target in order to protect 
it from deterioration. This was done by sealing the end cap 
which presses the target disc directly against the 0-ring 
vacuum seal, with a plate into which short lengths of tubing 
had been brazed, (See Figure 3.5) V/ater’ forced through 
the cap removed the heat from the target disc. 

3.3 The Neutron Source 



As in previous experimental measurements of neutron 
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Figure 3.5 Target Cooling System 
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slo'.ving do’.vn time (1,2), the Li*^(p ,n)Be'^ reaction v/ac used 
as the neutron source mainly on the "basis of its high 
yield and relatively good homogeneity. For a given, angle 
from the target and for laboratory proton energy up to 2.3^ 

MeV, the reaction is monoenergetic . The threshold of the 

7 7*^ 

Li'(p,n)Be' reaction is shov/n in Figure 4.6 of Adamantiades ’ 
thesis. Because of neutron intensity considerations,' a 
proton energy of. 2. 4 HeV was used in this experiment. 

The target used in these experiments consisted of a 
layer of natural lithium 0.8/cm^ thick deposited through 
evaporation on a circular tantalum disc 0.055 inches thick 
and 13/16 inch in diameter. Although the thickness of the 
target is less than 1/10 of the range of 2.5 MeV protons 
in lithium, some protons will undergo slowing down colli- 
sions before undergoing nuclear reactions. The energy 
spread of the actual beam, of protons incident on the tar- 
get, however, is very small in the Rockefeller Van de Graaff 
accelerator. 

To summ.arize, there were three major sources of neutron 
energy inhomogeneity in this experiment: 

1) . Angular dependence of neutron energy. 

2) Excited product nucleus reaction. 

3) Energy spread of the reacting protons. 

As mentioned in section 2.1, this energy inhomogeneity 
could be tolerated since the slowing down time is relatively 
insensitive to the source energy. 
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The i’.loderat injj; Assembly 

The moderating assembly (see Figure 5*'^) consisted of 
reactor-grade graphite stringers, approximately 4.0 inches 
by 4.25 inches and 4.0 feet long, stacked so that alternate 
flayers were at right angles to each other. The overall 

J 

dimensions of the assembly was 3.6 feet by 4.0 feet by 

I 

4.0 feet. The entire assembly was covered with aluminum 
foil to shield the monitoring and detection systems from 
R.F. pickup coming from the pulsing system, although this 
pickup was never entirely eliminated. 

In a stringer located at the center of one face of 
the pile, a hole was bored so that the target, along with the 
cooling water tubes and the cables of the monitoring sys- 
tem, could be inserted into the center of the assembly. 

'At a right angle to, and on the same horizontal plane 
as the beam tube was an open channel, 4.0 inches wide, 

4.25 inches high, and 19 .3 inches long, which led from the 
target to one face of the assembly. It was through this 
channel that the detector, foil, and graphite spacers were 
inserted or removed (see Section 4.2). 

3.5 Electronic Systems 
3 . 3.1 Beam Pulsing System 

Although the Rockefeller accelerator has a terminal 
pulsing capability, its single sweeping frequency was much 
too high for this experiment. Therefore, a post-accelera- 
tion system was used which had the advantage of a variable 
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pulse width and repetition, rate (see Figure 5*5)* A 
detailed discussion of the system is given in Reference (l). 
A repetition rate of 2kc/sec was used, allowing 500 usee 
between pulses, which was enough time to allow all traoi- 
sient neutron levels to die out. A pulse width of 200 nsec 
was used in t.his experiment. 

The pulsing system was triggered by the ^General Radio 

J 

wave generator (type no’. 1217-A) , which was also used to 
start the timing circuits of the analyzing system (see 
Section . 5 • 5 • 3) • The output pulse of the G.R. pulser' 
triggered an 3-H Research Laboratories Model No. 130 pulser. 

The variable width output pulse of the E-H pulser was 
the input for the driver circuit v/hich is essentially a 
pulse amplifier. An extensive discussion of the driver 
circuit is given by A.E. V/altar’ (3)> who designed it. ■ 

'V’aughan modified the circuit somewhat to give a narrov/er 
pulse. In this experiment the driver circuit delivered a 
2.9 ReV pulse to the deflection plates. 

3.5.2 The Detecti.on S:/stem 

I 

■ The object of these experiments v/as to detect neutron 
capture gamma rays as a function of time, so an efficient 
gamma-ray detector was required. A Harshaw Integral Line 
Scintillator Assembly, consisting of a 3 inch x 3 inch 
Nal(Tl) crystal attached to a 6363 Dulvlont photomutiplier 
tube, was used. The crystal and tube are hermetrically 
sealed within an aluminum housing, making the detector 
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a siini. le plug-iii A thermal neutron shield, consisting 

of 3 nim, thick boron carbide enclosed by aluminum, v/as pro- 
vided to' keep activation of the sodium at a minimum. 

The solid state preamplifier, v;hich plugged into the 
^ end of the detector assembly, was designed and built in a 
cylindrical can having the same diameter as the detector. 
Vaughan gives a schematic diagram of the preamplifier in 
Figure 4.6 of his thesis. 

The high voltage for the phototube v;as provided by 
a 2 KV power supply, while the preamplifier v;as powered by 
a small solid state power supply. The entire detector 
'was wrapped in -aluminum foil, and all leads to and from 
the detector were enclosed within copper braid to eliminate 

ii 

R.P. pickup. 

I 

5 . 3.3 The Analyzing; System ' • • * 

The main features of the time analyzing system we^re: 

1 ) A start pulse providing a time reference, derived 
from the beam pulsing- network. 

2) A stop pulse indicating the time at which an event 
of interest took place, derived from the detector. 

3) A timing unit, used in conjunction v/ith a multi- 
channel analyzer, which produced a count in the channel 
corresponding to the time interval betv/een the arrival of 
the start and stop pulses. Conventional, commercially 
available equipment v;as used in the entire analyzing system 
(see Figure 3*6). 
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Figure 3«o Analyzing System 
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The start pulse was obtained from the General Radio 
trigger (see Section 3.5-1). This pulse was fed through 
a General Radio 1392-A Time Delay Generator which provided 
a stable, variable time delay. Unf ortunately , the .output 
pulse of the G.R. delay was too narrow to trigger the tim- 
ing unit^ and an integrating circuit was employed to 
stretch the pulse. 

The detector output signal v/as fed to two Ortec 410 
Strobbed Discriminator, with the strobe input coming from 
an Ortec 407 Zero Cross-over Pickoff Unit. The cross- 
over pickoff was fed by the other Ortec 410, whose time con- 
stants had been optimized for a stable cross-over point. 

I 

I 

A positive logic signal was provided by the discriminator ! 

ii ' 

'when the incoming detector signal met minimum energy require- 

I 

ments. The strobe assured a constant time relation between 

I 

the detector input and discriminator output. The positive | 

I 

logic pulse was attenuated by a resistor netw'ork before it j 
was fed into the timing unit. I 

I 

The timing unit consisted of a TMG Model 211 Time of | 
Plight Logic Unit v/hich plugged into the TMG Model CN-110 ' 

j 

256 Channel Pulse Analyzer. The 211 unit converted the time ; 
difference betv/een the trigger input and a detector signal i 

I 

to an approoriate channel number, depending on the scale I 

i 

factor which could be varied from 0.25 to 6-4. psec/channel . i 
This information v;as then stored in the multi-channel analy- , 
zer. ' ! 

The TMG 211 Unit could handle more than one stop pulse , 

i 

( 
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(detecoor signal) per cycle. However, 16 psec v;ere 
needed to analyze ainl record each incoming signal. A 
correction for this deadtime was included in the computer 
program used to process the data (see Section 4.3). More 
detailed information on the 211 unit and the multi-channel 
an,alyzer is given in Reference (4). 

3.3.4 The ?;lonitoring and Normalizing System 

To insure that the pulsing system v/as operating under 
the best possible conditions at all times a beam pulse 
monitor v/as provided. 

The beam pulse monitor v/as necessary due to the pro- 

I 

blem of focusing the proton beam from the accelerator. ' 

Because of various changes and drifts in the accelerator, • 

the setting for optimum focusing changed from day to day, i 
and had to be constantly adjusted. In pulsed operation, 

I 

the only practical method of focusing was to use an electri- 1 

I 

cal signal from the target itself and to observe the signal ! 

! 

with a fast oscilloscooe . The system employed is sketched 

I 

in Figure 3*7. Notice that the target was insulated from : 

the beam tube with a glass ring which enables the signal ^ 

to be picked off the target. The monitoring system was | 

I 

esrecially susceotible to R.F. pick-up from the pulser. ' 

* "• • i 

All cables of the system were shielded against pick-up j ‘ 

but the 'problem was never completely eliminated. A typical 

j 

target pulse can be seen in Figure 4.11 of Vaughan’s thesis. • 

I 

The electronic pulse being delivered to the deflection plates^ 
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Pigure 5*7 Monitoring Systems 
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could al.'30 08 Eionitorsd. This systeia is also shov/n in 
Fi^re 5-?. 

A neutron monitor was used as a means of normalization 

for comparing different runs and subtracting background. 

A detailed discussion of normalization is given in Chapter 

5 of Reference ( 1 ). As in previous experiments, this author 

concluded that total neutrons produced during an experiment 

would be the most reliable quantity to use of normalization. 

A boron triflouride counter was placed between the graphite. 

stringers of the assembly about 13 inches from the target 

and on the side opposite the absorber foil and detector. 

Thus, with the output ‘of the BP^ tube going to a scaler, 

the time integrated neutron -flux as seen by the detector 

t) 

could be recordci. 

I 

i 
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Experimental Results and Discussion 

A-. 1 T he Use of Resonance Absorbers as Neutron Density 

, Indicators 

i 

Resonance absorbers were used in this experiment to 
sample neutron population- at the resonance energy. The 
most suitable materials for this purpose are those which 
exhibit a strong isolated absorption resonance. The 
higher and narrower the resonance, the better it approxi- 
mates a delta -function. Table II of Reference (1) gives 
a list of absorbers which exhibit strong isolated reson- 
ances. The energy and peak cross section of the reson- 
ance' is also given. 

Another important consideration in choosing a suit- • 
able absorber was its prompt neutron capture gamma-ray 
spectrum. This did not prove to be a problem since the 
de-excitation gamma rays of most materials are spread 
over a wide range of energy. Reference (1) gives capture 
gamma-ray spectra for the three absorbers used in these 
experiments: cadmium, indium, and gold. The capture gamma 

spectra of these materials extends from the KeV region ' 
to energies exceeding 7 "v-eV. = 

Indium and gold were chosen for their strong isolated, 
resonances at l.h6 eV and 5*0 respectively. Cadmium, 



.vhich has a hroad I’asonaace that reaches a maximum at 
about 0.2 eV, was also used, although it v/as expected that 
the time resolution for these measurements v/ould not be 
as good as that for indium and gold. 

'd.2 Performance of the Experiment 

In order to insure a constant geometric relationship 
between the foil and detector, a foil holder was used. 

The foi'l holder consisted of a one inch thick slab of gra- 
phite, which v/as^sized to fit into' the channel of the assem 
bly and was bonded to the end of the detector with epoxy 
cement. The foil v/as taped to this slab, and the v/hole 
detector-foil array was inserted into the graphite assem- 
bly thi'ough the channel (see Figure -T.l). u 

■ To measure the spatial dependence of the slowing down 
time, it was necessary to be able to move the detector- 
foil array along the channel. Graphite slabs of varying 
thickness, which v{ere inserted in the channel between the 
source and foil, assured that there would always be solid 
moderating material in front of the foil. 

Before performing the actual, experiments, short- runs ■ 
were made while varying the discriminator of the Ortec 
415 Strobbed Discriminator. For each foil, the discrimin- 
ator setting which resulted in the best foil counts to 
background .counts ratio v/as determined. 

A typical set of experimental data is shov/n in Figure 

N 

> 

4.2(a). .These results are for a 7 cm x 7 cm x 0.05 cm gold 
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Figure 4ol Horizontal Mid-^Plane of Moderator Assembly 



foil pooitibood 23 cm av;ay from the target. The beam 
pulse •.vidL.h was 200 nsec,. and the time interval betv/een 
bursts ’.vas 500 usec. The channel width was 2^0 nsec. 

A time distribution with the gold absorber in place was 
. ta.ken, and a distribution without the absorber v/as taken 
immediaP^elj aft'erwards . The BF, counts for the gold foil 
was 251,111 and for the background was 259,992. Notice 
the very prominent peak around channel 20. This peak is 
called t.he prompt peak and occurs when the burst of pro- 
tons hits the lithium target. This peak was used as an 
''•accurate zero time reference in the experiment. The right 
half of the graph (notice the change of scale) shows the 
region in which the most probable slov/ing down time occurred. 

t! 

No peak is discernible from the raw data. 



4.3 Normalization, Backa;round Subtraction, Dead Time 

Correction and Smoothing 

As mentioned in Section 3* 5*^, the time integrated 
neutron flux, as derived from a BF^ tube, was used for 
normalization. The number of BF^ counts registered by 
the scaler v/as recorded for each experiment. A factor 
obtained by taking the ratio of background-run BF^ counts 
to foil-run BF^ counts was applied to the foil data. 

In order to obtain the net effect of neutron cap- 
tures in the foil, the counts originating from all other 
events must be subtracted when the time distribution is 
taken. These sources of parasitic gamma rays include: 
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Figure ^.2 (a) Experimental Data From a Background 
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1) Inelastic scattering by the target, the absorber 
foil, and the detec lox'. 

2) • Absorption in materials other than the absorber 

ft 

foil. 

5) Reactions coincident with the neutron burst. 

I 

4) Overlapping events from preceding bursts. 

5) Steady state background. 

All these items were taken into account by taking 
a time distribution without the absorber foil in place, 
and then subtracting it from the distribution taken with 
the foil. 

As mentioned in Section there was a 16 psec 

dead time associated with each count registered 'in the 
multi-channel analyzer. In Appendix C of Reference (2), 
a theoretical discussion of the dead time correction for 
the Ti'.'C Model 211 Unit is presented. The correction fac- 
tors derived in that reference were used in this thesis. 

Computer programs v/hich normalized, subtracted the 
background, and applied dead time corrections to the raw 
data are listed in Appendix A. 

Figure 4.2(b) shoves the time distribution resulting 
from processing the raw data in Figure 4.2(a)". It is 
obvious from this graph that there is a maocimum in the 
time distribution occurring roughly at channel number 60. 
It is also obvious, however, that the fluctuations in the 
curve prevent an accurate determination of the maximum 
point. To reduce this uncertainty, a computer program 
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Figure ^.2 (b) Time Distribution . at 5*0 eV and 28.00 cm Away From the Source 
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developed at iVI.I.-. (3) v;hich applied a s.^ioothing tech- 
nique to the corrected data was used. Effectively, the 
program -fits a Fouier series curve to the data points. 
Smoothing is accomplished by choosing not to fit the ser- 
ies tojir.igh frequency waves, which correspond here to sta- . 

iP\' 

tistioMi fluctuations. The frequency . cutoff , above which 
all h:-|her frequencies will be thrown away, is under the 
contri'^’ of the programrasr. Thus, any desired degree of 
smoothing can be effected to a set of data points. In 
this experiment, all frequencies less than four channels 
wide were dropped. Figure 4.2(c) is the smoothed version 
of Figure 4.2(b). Using this smoothed distribution j it 
was estimated that the maximum point of the time dis- 

ii 

tribution occurred in channel 62.5+1. Since the prompt 
peak is in channel 20, the slowing down time is: (43.5j;l) 

X (.25 psec/channel) = 10.8?+. 25 psec. The smoothing 
program is listed in Appendix A. 



4.4 Space and Energy Uenendent Slowing Down Time Mea- 
surem e nts 

If various absorbers with resonances located at 
different energies are used and if the distance of these 
absorbers from the source is varied, the time dependent 
neutron distribution can be measured as a function of 
space and neutron energy. 

Figures 4.3 to 4.6 represent time distributions to 
the cadmium resonance^ Figures 4.7 to 4.11 represent time 
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CHRNNEl. NUMBER 

Figure 4.2 (c) Smoothed Version of Figure ,4.2 (b) 



dint ion. to the indium reoonance' and Figures 4.15 
to 4.16 rexiresent time distributions to the gold reson- 
ance. In each Figure, (a) repx''esents the distribution 
before smoothing and (b) after smoothing. Notice that 
the scale for the cadmium curves is 500 nsec/channel , 
Vvhile it is 250 nsec/channel for the indium and gold 
curves. 

The-' results of this experiment, together with their 
probable errors, are graphically presented in. Figures 
4.17 to 4 . 19 . The experimental uncertainty associated 
with each point reflects the broadness in the time dis- 
tribution- around the maximum value. As expected, the 
peaks of the cadmium distribution were much broader than 
those of gold or indium'. The solid line in the gold and 
indium figures has been constructed from theoretical time 
distributions' (see Section 2.5), based on diffusion theor 
(4). 

The experimental results for. the slowing dovm time 
to 5*0 eV appear to be in good agreement with the theo- 
retical calculations. The agreement v/ith theory for 
1.46 eV is not as good. The experimental results indi- 
cate that there is little, if any, spatial dependence of 
the slowing down time for distances closer than about 
50 . cm from the source. For distances greater than 50 cm, 
the ex-oerimental results show a sharp rise in t^ with dis 
tance, indicating a strong spatial dependence. Theore- 
tical calculations do not predict as sharp a rise of 
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Figure 4.5 (a) Time Distribution at 0.2 eV and 12,5 cm Away From the Source 
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Figure 4.4 (a) Time Distribution at 0.2 eV and 18.0 cm Away From the Source 
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Figure 4.4 (b) Smoothed Version of Figure 4.4 (a) 
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Figure 4.5 (a) Time Distribution at 0.2 eV and 28.0 cm Away From the Source 
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Figure 4.5 (b) Smoothed Version of Figure 
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Figure 4.6 (a) Time Distribution at 0*2 eV and 44.5 cm Away From the Source 
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Figure 4-. 7 (a) Time Distribution at 1.46 eV and 12.5 cm Away From the Source 
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Figure 4,7 (b) Smoothed Version of Figure 4.7 (a) 
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Figure 4,8 (a) Time Distribution at 1.46 eV and 18.0 cm Away From the Source 
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Figure 4.8 (b) Smoothed Version of Figure 4.8 (a) 
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Figure 4.9 (a) Time Distribution at 1.46 eV and 28.0 cm Away From the Source 
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Figure 4.9 (b) Smoothed Version of Figure 
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Figure 4.10 (a) Time Distribution at 1.46 eV and 57*5 cm Av/ay From the Source 
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Figure 4.10 (b) Smoothed Version of Figure 4.10 (a) 
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Figure 4,11 (a) Time Distribution at 1.46 eV and 44.5 cm Away From the Source 
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Figure 4.11 (b) Smoothed Version of Figur 
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Figure 4.13 (a) Time Distribution at 5*0 eV and 12,3 cm Away From the Source 
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Figure 4.15 (b) Smoothed Version of Figure 4.13 (a) 
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Figure 4.14 (a) Time Distribution at 5»0 eV and 18.0 cm Away From the Source 
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Figure 4.15 (t>) ' Smoothed Version of Figure 4.15 (a) 







1-igure -4-. 16 (a) Time Distribution at 5.0 eV and 44.0 cm Away From the Source 



as sho'.vii bj' the exueriraents . 

■ This .disagreeneiifc between theory and experiment v/as 
probably due to the breakdov/n of the infinite medium 
bo\mdary conditions, since the foil was only 9 cm from 
the outer edge of the assembly at the last spatial point. 
The absence of moderating material surrounding the foil 
meant that the neutrons took a longer time to slow down 
to a given energy. More measurements need to be taken 
at greater distances to resolve this question. 

A few measurements were made with the detector separ- 
ated from the foil ’ 03 ;- about 10 cm of graphite in order 
to determine if the perturbation caused by the detector 
had any effect on the slov/ing down time. These '.results 

% ■ 

are marked by triangles in the Figures. It can ii>e seen 
that the gold and cadmium results agree, v/ithin the limits 
of experimental uncertaintity , with previous measurements. 
The indium point, however, falls far- below the original 
measurements. The author can give no explanation for the 
contradiction of this point with either the original result 
for indium or the close agreement betv/een the first and 
second measurements for the gold and cadmium. Further 
measurements, with the detector and foil separated, need 
to be performed to clear up this contradiction, 

h .3 Concltisions 

There was some agreement between the spatially 'depen- 
dent slowing down time measurements made in these experi- 
ments and theoretical calculations using diffusion theory. 
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Figure 4.17 Slowing Down Time vs. Distance at 



Figure ‘4-, 18 Slowing Down Time vs. Distance, at 1.-4-6 eV 
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Figure 4.19 Slowing Down Time vs. Distance at 5.0 eV 
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No measurable spatia.1 dependence of the slowing down time 
was observed for distances less than 50 cm. Spatial depen- 
dence was observed for distances greater than 50 cm, although 
this depend‘ence was much more promounced in the experiments 
than predicted hj slowing dovm theory. More space depen- 
dent measurements at distances greater than cm should 
be made in a larger graphite assembly to determine the 
true spatial dependence. Measurements to determine ’if 
the perturbation of the detector had any effect on the 
slowing dov/n time were inconclusive and should be repeated. 
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Appeniix A. Coinputor Program Listings 



MORMALI ZATI0\S DEAD TIME CORRECTION, AND BACKGROUND 
SUBTRACTION CODE LISTING 

DIMENSION TBM( 2 55 } , CHAM ( 255) , B ACK { 255 ) , FO I L I 25 5 ) , ANORM ( 255 ) , 
1CCACK(2351 ,CFOiL(255) ,X(255) , Y ( 255 ) , T8A ( 2D6 ) , ACHAN ( 206 ) ,B{800) 
COMMON T8M, CHAN, WIDTH 

EQUIVALENCE ( IBM ( 50 ) , T BA ( I ) 1 , ( CHAN C 50) , ACHAN 1 1 ) ) 

READ 200, WIDTH 
200 F0RMAT(F6.2) 

READ 99,aC0UNT,FC0UNr 

99 FORMAT (2F10.0) / 

NL!M = 0 . ' / 

103 NUM=NUM+1 

READ 100, (8ACK( I ), 1=1,255; 

100 FORMAT (3X,10F7.0) 

READ 101, (FOIL ( I ), 1=1,255) 

101 FORMAT (3X,i0F7.0) 

READ 98,3BURtABUR i 

98 FORMAT (2F10-0) j 

PRINT 701,8COUNT,FCOUNT i 

701 FORMAT (3X,2F10.0) 

PRINT 100, (BACK( I ) ,1=1 ,255) 

PRINT 101, (FOILIl ) ,1=1,255) ! 

PRINT 702,86UR,ABUR 

702 FORMAT 13X,2F10.0) 

L = -l 1 

420 L = Lt-l 

IF (L)419,419,417 ; 

419 00 418 1=1,255 | 

418 X{ I )=BACK(I ) 

BUR=BBUR I 

GO TO 415 ' 

417 DO 416 1=1,255 
416 X( I )=FOIL(I ) 

6Ua=ABUR : 

. 415 0UM=L6. /WIDTH 

M=DUM . ; 

Y( l)=BUR*L0GF(BUR/{8UR-xm) ) | 

SUH=0. , 

DO 610 1=2, N I 

SUH=SUM+X( I-l ) • I 
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610 Y( I)=BUR9L0GF( ( BUR-SUM) /(BUR-SUM-XII ) ) ) 

SUM^SUM^Xi.Nl 
NM = .'^+L 

‘ ■ DQ 615 l=MM,255 

Y( I } =BUr<RL0GF( (BUR-SUM) /IBUR-SUM-X(I )) ) 

, K=I-M 

615 SUM=SUM+X( I )-XtK) 

IF { L ) 414, ^12 

414 DO 413 1=1,255 
413 CBACKt I )=Y( I ) 

GO TO 420 

412 DO 411 1=1,255 
411 CFOIL( I )=Y( I ) 

2 = BCOUNT/FCOJ.NT 
DO 906 1=1,255 
ANORM( I ) =Z»CFOIL{ I ) 

TBMd ) =ANORM( I 1 -CBACK( I ) 

CHA:M( [ )=FLOA‘rFl I ) 

906 CONTINUE 

PRINT 306, I T6M( I ) , 1 = 1,255) I, 

306 FORMAT IlH ,10F10.1) 

CALCOMP ROUTINES 

50 CALL PLOTSl (B,800) 

CALL SYK3L5(3.0,9.0,0.2,20H ID. NO. = M3994-5788 , 0 . , 20 ) 
CALL N1JM8RU4.0,8.5,D.2,NUM,0, ,-U 
CALL P ICrURI 8.5,4.2, 15H CHANNEL NUMBER, 

115, 19H COUNTS PER CH ANNEL , 19 , ACHAN, TBA , 206 , 0. , 0 ) 

END FILE 10 

READ 99,3C0UNT,FC0UNT 
IF { BCOUNT) 700,700, 103 
700 CALL EXIT 
END 
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data 5M0JTHING CODE LISTING INCLUDING NORMALIZATION, DEAD 
TIME CORRECTION, AND BACKGROUND SUBTRACTION 

DIMENSION TBN( 37) ,XSINI75) ,SEIGEN(37,75) ,XC0S1 (75),XCOS2(75) , 
1XC0S4(37 J , IBM (255) , TDK (255 5 , CHAN (255) , THE! A ( 75 5 , ALPHA ( 75 ) , 
2BETA{75) , GAMMA (75 5 ,WT(75) , B ( 800 5 , BACK ( 255 ) , FOI L ( 255 ) ,ANDRM(255), 
3CBACK(255) ,CFOIL(255),X(2555 , Y ( 255 ) , TBA ( 21 6 ) , ACHAN ( 2 1 6 5 
COMMON TBN,XSI>I,SEIGEN,XC0S1,XC0S2,XCDS4,TBM,T6K, CHAN, 'WIDTH 
EQUI VALENCE (TBK{ 40) , T6A( 1) } , ( CHAN ( 40 5 , ACHAN ( 1 ) ) 

READ 200, IH,OMEGAH,CRT,FACTR, WIDTH 
ZOO FORMAT ( I 4, E 10 . 4, F6. 0 , F6. 3 , F6. 2 ) . • 

READ 102, ( WT( I ) , 1=1, IW) 

102 FORMAT ( 12F6.3) 

PRINT 201, IW,OMEGAH, CRT, FACTR, WIDTH 
201 FORMAT (1H1,9H N-OMEGA= , 1 4 , IX, 9H D-OMEG A = , E 10 . 4, 1 X , 

15H CRT=,F6.0, 1X,8H F ACTOR= , F6. 3 , 7H WIDTH=,F6.2) 

PRINT 205 

205 FORMAT ( IH ,20X,16H WEIGHT FUNCTION) 

PRINT 203, (WT( I ) , 1=1, IW) 

203 FORMAT (IH ,5615.7) 

READ 99,BC0UNT,FC0UNT 
99 FORMAT { 2F10.0) 

MUM=0 / 

103 NUM=NUM+1. • ■ • ‘ 

^ READ 100, (BACK(I) , 1=1, 255) 

100 FORMAT (3X,10F7.0) 

READ 101, (FOIL! I ), 1=1,255) 

101 FORMAT (3X,10F7.05 
. READ 98,BBUR,ABUR 

98 FORMAT {2F10.0) 

PRINT 701,3C0UNT, FCOUNT 

701 FORMAT (3X,2F10.0) 

PRINT 100, (BACK! I ), 1=1,255) 

PRINT 101 , (FOIL( I ), 1=1,255) 

PRINT 702,B8UR, ABUR 

702 FORMAT (3X,2F10.0) 

L=-l ' . 

420 L=L+1 

IF IL)419,419,417 
419 DO 418 1=1,255 
418 X{I)=8ACK(I) 

BUR = BB'JR 
GO TO 415 

417 DO 416 1=1,255 
416 X( I )=FOIL( I ) 

BUR=ABUa 



I 



I 



I 
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0U.M = I6. /WIDTH 
N^DUM 

YC1)=BUR"L0GFIBUR/(3UR-X(l) ) ) 

SUM=0. 

DO 610 1=2, N 

SLIM = SUM + X( I- I) 

610 Y( i )=3UR=L0GFC (BUR-SUM)/{QUR-SUM-X(I ) ) ) 

SUM=SUH^X(N5 

NM=N<-1 

‘ 00 615 I=NM,255 

• Y( I )=BUR*LOGF{ { BUR-SUM )/( BUR-SUM-X II ) ) ) 

K=I-N ■ / 

615 SUM = SUM-»X( n-X(K) 

IF (1)414,^14,412 
AL4 on 413 1=1,255 
413 CBACK{ I )=Y( I ) 

GO TO 420 

412 DO 411 1=1,255 

411 CFOIU I )=Y( I) I 

DO 905 1 = 1,37 ; 

YKAKU4=0M6GAH*FL0ATFIIW-1) i^FLOATFI I ) • ■ 

XC0S4( I 1=C0SF(YKAKU4) 

905 CONTINUE 

DO 901 1=1,37 ' 

00 902 J=1,IW 

XKAKU=OnEGA.H^FLOATF ( I ) »FLOATF( J-1 ) : 

SEIGENI T , J 1=SINF(XKAKU) t 

902 CONTINUE | 

901 CONTINUE | 

DO 903 j=i, iw ; 

THETAI J i =0MEGAH»FL0ATF I J-1 ) ' 

ALPHA{ J ) =THETA{ J ) =•* { -3 ) * ( THETA ( J ) *»2 + THETA ( J ) *SE IGEN (1 , J ) • 
1C0SF(THETA( J ) ) -2. »SE I GEN ( 1 , J ) * *2 ) i 

BETA( J) = 2,»THETA( J)»^»(-3)»(THETA( J)*(l,+COSFtTHETA( J) )«*2)-2.* 
1S£IGEN( 1,J ) »C0SF{IHETA£ J) ) ) 

GAMMA! J)=4.i TH£rAlJ) »i^{-3)»{ SEIGEN( 1,J)-THETA(J) »COSFtTHETA(J) ) ) 

903 CONTINUE 

DO 906 1=1,255 

2=SC0UNT/FC0UNT ' , j 

ANORMI I)=Z-’^CFOIU I) 

TBM( I)=ANORM( I )-CBACKm ' ! 

906 CONTINUE ' ' 

DO 15 1=1,27 ! 



I 



I 
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00 90't J = 1,IK 
YKAKUl=nMEGAH»FLQATF ( J-1 ) 

'• ■ Y!<AKU2 = 27. *YKAKUl 
' XCOSl ( J ) =C0SF { YKAKUl ) 

XC0S2{ J)=COSF( YKAKU2) 

90 A CONTINUE 

DO 511 1=1, I W 

.4;„ YA = ALPHA(l) . ; 

YB=OCTA(I) 

YC = GAM<MAl I ) 

SIG4=0. 

SIG3=0,5* ITBM{1)*S£IGEN{U n+TBN(27)«SEIGEN{27,I ) ) 

DO 1 NN=3,25,2 

1 SIG3 = SIG3i-TBN(NN) sSElGENINN, I) 

00 2 NN=2,26,2 

2 SIGA=SIGA+TBN(NN)*SEIGEN(NN, I> ' 

XA = YA<^(T8N(1)»XC0S1( I ) -TBN t 27 ) »XC0S2 ( I ) H-YB*S I G3 + YC»S I GA 
xsrm )=XA*WT( r ) 

511 CONTINUE II 

DO 517 IX=1,16 
YA = ALPHA( IXt-1) 

YB = BErAHX^-l ) 

YC = GAHMA{ IX+U 
SIGA=0. 

SIG3=0.5«ftXSIN( 1) *S£IGEN( IX, U + XSIN( I W ) *SE I GEN ( IX, I W ) ) 

NY=I W-2 

DO 3 NN=3,NY,2 

3 SIG3 = SIG3 + XSIM(NN)«-SEIGENHX,NN) 

NY=IH-i 

DO V NN=2,NY,2 

4 SIG4=SIG4+XSIN(NN)»SEIGENI IX,NN) 

XA = OMEGAH» ( YA* ( XS I N ( I ) -X S I N { I W ) * XC0S4 ( I X ) ) YB» S I G3 + YC*S I G4 ) 
IN=IX 

T8K1 IN)=0.63662*XA 
CHA'UIM) =FLCATFI IN) 

“ 517 CONTINUE 

DO 914 J=1,IW 

YKAKUl = OKEGAH*FLOATFU-l) 

yKAXU2=37,»YKAXUl 

XCOSK J)=C0SF(YKAXU1) 

XCOS2( J)=C0SF(YKAKU2) 
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914 CONrnUE 

DO 530 J=2,15 
DO 525 1=1,37 
K=16*( J-l)-10<-I 
525 TSN(I)=TBM(K) 

DO 521 I = U IW 
YA=ALPHA(I) 

YB=BETAII) 

YC=GAMMA1I) 

SIG4=0. 

5 103 = 0.5* (TB.'H 1 ) ^ SE I GEiN: ( 1 , I ) +TBN { 37 ) *SE I GEN { 37 , I ) ) 

DO 5 NN=3,35v2 

5 SIG3=SIG3+TBN(NN)*SEIGEN(NN,I) 

DO 6 NN=2,36r2 

6 S1G4=SIG4+TBN(NN)*S£IGEN(NN, I ) 

XA = YA«f [TBN{ 1 ) sXCOSU I ) - T BN ( 37 ) *XC0S2 ( I ) ) +YB»S I G3 + YC»S I G4 
xsiN{ I )=xA*wrn ) 

521 CONTINUE 

DC 527 IX=11,26 
YA=ALPHA( IX + 1 i 
YB=BETA( IX+L) 

YC = GAMMA{ IX + U 
SIG4=0. 

SIG3=0.5>(XSIN{ 1)»SEIGENIIX,IH-XSINI IW)»SEIGEN{ IX, IW) ) 
NY=IW-2 • 

DO 7 NN=3,NY,2 

7 SIG3 = SIG3<-XSIN(NN)*SEIG£N(IX,NN) . 

NY=[W-1 

DO 8 NN=2,NY,2 

8 SIG4 = S.IG4 + XSIN(NN)*SEIGEN(IX,NN) 

XA = DMEGAH* ( Y A « ( XS I N ( 1 ) -XS I N I I W ) »XC0S4 ( IX J ) +YB*S I G3<-YC*S I G4 ) 
IN=16» U-1 )+IX-10 
TBXl IN)=0.63662»XA 
CHANi IN) =FLOATF ( IN) 

'527 ' CONTINUE 
530 CONTINUE 

DO 535 1=1,25 
K=Z30+I 

535 T8NC I ) =TBM{ K) 

DO 924 J=1,IW 
YKAKU1=0MEGAH»FL0ATF(J-1) 

YXAKU2 = 25. »Y:<AFJJi 
XC0S1{J)'=C0SF( YKAKUl) 

XC0S2I J)=C0SF{YKAKU2) 
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92 CCNTINUC- 

DQ 531 I=i,IW 
YA=ALPHA( I } 

YB=B£I A( I ) 

YC = GAMfSA{I) 

SIC-'» = 0. 

SIG3 = 0.5*( rONl 1 ) «S£IGENt 1, I ) +T BN ( 25 ) »SE I GEN { 25 » I ) ) 

DO 9 NN-3,23,2 

9 SrGi=SIG3^ TBN{NN)ftS£iGEN(NN, n 
DO 10 MN=2,2A,2 

10 SIG4 = SIG4 !-T3N(NN5 »SEIGEN(NN. I ) 

XA = YA» ( TBN { 1 ) «XC0S1 ( I ) -TBN{25) «-XC0S2 ( I ) ) +YB*S I G3 + Y S I G4 
XSI'H I )=XA»WT{ I) 

531 CONIINUE 

DO 537 IX=11,25 
YA^ALPHA'i IX+1 ) 

YB = B£TA( IX + 1 ) 

YC=GAMHA(IX+1) \ 

siG4=o. ; 

SiG3 = 0.5K tXSIN(U«-SEIGENl*IX,l)+XSIM( IW)«^SEIGEN( IX, IW) ) 
NY=IW-2 

DO U NN=3,MY,2 

11 SIG3 = 5!G3+-XSIN(NN)^^SEIGEN(IX,NN5 
MY=IW-1 

00 12 MN=2tNY»2 1 

12 SIG4-SIG4-t-XSIN(NN)*SEIGEN( IX,NN) i 

XA=OMEGAHfr lYAitlXS!N{ 1 )-XSINC I W } *XC0S4 ( I X ) 1 +YB*S I G3 + YC*S I G4 ) 
IM=230+IX [ 

TBK( IN)=0-63662»XA 

CHAN{IN)=FLOATF( IN) 1 

537 CONTINUE . 

PRINT 307 ! 

307 FORMAT {1H1,18H SMOOTHED FUNCTION) 

PRINT 306, ( TBxn } , I = l»255) 

306 FORMAT I IH tIOFIO.I) ! 

CALCOMP ROUTINES 

• ■ I 

50 CALL PLOTSKB'jBOO) 

CALL NUMBRH4,0,8.5,0.2»NUM,0. ,-l) 

CALL SYf'.BLSI 3.0,9,0,0.2,20H ID NO. = H3994-5788 , 0 . , 20 ) 

CALL PICrUR{6.5,4.2, 15H CHANNEL NUMBER, ! 

7 115. 28H SMOOTHED COUNTS PER CHANNEL , 28 , ACH AN , T B A , 2 1 6, 0 . , 0 ) 

END FILE 10 , i 

c : 

READ 99,BC0UNT.FCCUNT ' i' 

IF laCOUNT) 700,700, 103 



700 CALL EXIT 
END 



I 
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